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The key requirement for harnessing the quantum properties of light is the capability to detect 
and count individual photons. Of particular interest are photon-number-resolving detectors, which 
allow one to determine whether a state of light is classical or genuinely quantum. Existing schemes 
for addressing this challenge rely on a proportional conversion of photons to electrons. As such, 
they are capable of correctly characterizing small photon fluxes, yet are limited by uncertainties 
in the conversion rate. In this work, we employ a divide-and-conquer approach to overcome these 
limitations by transforming the incident fields into uniform distributions that readily lend themselves 
for characterization by standard on-off detectors. Since the exact statistics of the light stream are 
obtained from the click statistics, our technique is freely scalable to accommodate - in principle - 
arbitrarily large photon fluxes. Our experiments pave the way towards genuine integrated photon- 
number-resolving detectors for advanced on-chip photonic quantum networks. 


Quantum information science is at the cutting edge of 
modern physics and technology. In this context, perhaps 
the most ambitious goal is to realize scalable quantum 
information processing and computing based exclusively 
on linear optical configurations and photon-counting 
devices CH3|. Notably, any such optical quantum¬ 
computing scheme hinges on the ability to detect and 
manipulate the states of light at the single-photon level: 
Quantum cryptography, entanglement swapping, and 
quantum teleportation, to name a few, would clearly 
be impossible without reliable single-photon-counting de¬ 
vices min]- Moreover, exact photon counts provide ac¬ 
cess to genuine photon number statistics, and in turn 
are the principal means of reliably establishing the non- 
classicality of any type of light field CHE]- Another 
potential application of photon-number-resolving detec¬ 
tors (PNRs) was recently highlighted in the context of 
coherent optical communications m, where they enable 
coherent optical communications with a performance su¬ 
perior to the standard quantum limit, even in the high 
mean photon number regime. 

Existing schemes for measurements at the single¬ 
photon level employ on-off detectors, e.g. avalanche pho¬ 
todiodes (APDs) [IT] . and as such are inherently limited 
by the so-called dead time. When an APD is triggered, it 
typically remains “blind” for several nanoseconds there¬ 
after, and as a result, succeeding photons impinging on 
the detector cannot be registered m- In addition to 
being detrimental to the overall detection efficiency, this 
effect may corrupt the very photon statistics one strives 
to characterize. Moreover, this saturation effect also in¬ 
troduces undesired correlations to the count sequences. 
To this date, the perhaps most noticeable scheme for 
PNRs is based on superconducting nanowires [20] . Yet, 
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on many occasions, cryogenic measurements may be im¬ 
practical, or the incident photon flux may exceed the 
capacity of the system. Evidently, a fundamentally dif¬ 
ferent approach will be required to reconcile the demands 
for high speed, low noise, and maximized quantum effi¬ 
ciency with the ever increasing count rates required by 
modern technologies [2TH25] . 

In contrast to PNRs, on-off detectors deliver well- 
defined “clicks” upon excitations with any non-zero num¬ 
ber of photons [26] . Consequently, the by far most acces¬ 
sible quantum-optical measurements are click-counting 
statistics, instead of actual photon counts m- The ques¬ 
tion naturally arises as to whether it is possible to cir¬ 
cumvent the limitations of on-off detectors, and to ex¬ 
ploit these robust and widely available components to 
accurately measure multiple photons. 

In this work, we propose, implement and characterize 
a photon counting device based on a multiplexed array 
of on-off detectors. In our arrangement, the discrete evo¬ 
lution dynamics of the incident light field is manipulated 
so as to distribute the photons uniformly between the in¬ 
dividual channels. Crucially, the click-counting statistics 
obtained from this type of multiplexed sensor converge to 
the actual photon-counting statistics as the number of in¬ 
dependent detectors is increased [28]. Along these lines, 
it becomes possible to reliably probe the non-classicality 
of arbitrary optical input fields in a readily scalable and 
integrated fashion. 

Let us first consider a stream of single-photon states 
being routed through a uniform l-to-TV multiplexer and 
onto an array of APDs (see Fig. [l] top). In our method, 
this is achieved by cascading m stages of 50/50 beam 
splitters, yielding N = 2 m . Under these premises, ev¬ 
ery single-photon will have a probability of 1/N to be 
detected in one of the N channels. Due to the spatially 
extended wave function, any two incoming photons are 
likely to be found in different outputs with a probabil¬ 
ity of 1 — 1/N: When one photon is detected, the global 
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FIG. 1: Top: Sketch of the l-to-8 optical integrated multi¬ 
plexer consisting of three beam splitting stages. The output 
fields are fed into APDs. Bottom: Output click statistic of the 
above multiplexer for classical attenuated laser light (blue), 
as well as single-photon Fock-states (green). In the experi¬ 
ment, both input states yield a flat uniform output statistic 
provided by the high quality of the optical-integrated device. 

probability of registering the next photon in any of the 
remaining APDs is (TV —1) times greater than in the same 
one. This remains true even if the two photons enter the 
system simultaneously. In this manner, the fidelity of the 
device improves with the number of output ports, and is 
even independent of the type of input state [28]. 

To experimentally demonstrate the functionality of our 
approach, we realized a discrete network of integrated 
50/50 beam splitters cascaded in m = 3 steps, yielding 
a total of TV = 8 output channels. These photonic struc¬ 
tures were implemented in fused silica glass by means of 
the femtosecond laser writing technique H EDI, see Ap¬ 
pendix [A] As input states we consider the two limiting 
cases: i) Low-intensity laser light, which can be idealized 
as a quantum coherent state [5]. As such, its statistics 
feature substantial yet temporally fluctuating bunching 
of photons, making it a perfect test case for the photon 
number resolving capabilities of our setup, ii) Heralded 
single photons from a spontaneous parametric down con¬ 
version source, which approximately represent the ideal 
scenario of single-photon Fock states [3T] , 

In the multiplexer, classical coherent states l^in) = |a) 
are split into eight spatially separated coherent states 
of equal amplitude, i.e. |Vwt) = \—a/y/8) |ia/\/8) 
\a/y/S) \ia/y/S) |— a/y/S) |i a/y/S) \—a/y/S) \—ia/y/S). 
Consequently, a perfect photon counting characterization 
should yield a Poissonian photon number distribution. 
When one instead evaluates the counting click coinci¬ 


dences, it can be analytically shown that the resulting 
click-counting statistics have to follow a binomial distri¬ 
bution [28] . Note that in case of a sub- (or super-) Poisso¬ 
nian photon number distribution, it likewise follows that 
the click statistics are sub- (or super-) binomial, respec¬ 
tively (see for instance Fig. [2]). 



FIG. 2: Theoretical click-counting statistics for three dif¬ 
ferent photon number distributions at an average number of 
clicks of 0.25 per measurement. The resulting values are 0 
for binomial, —0.22 for sub-binomial, and +0.22 for super¬ 
binomial click statistics. 

In general, this behavior with respect to TV on-off de¬ 
tectors is mathematically described by the expectation 

value Ck = (:(£K fe (i -7t) N ~ km ^ Here, the normal 

ordering for the bosonic operators [32] is indicated by 
the symbol the number of APDs that click in a 

certain time window is 0 < k < TV, and the operator 
7r — :1 — exp[—(ryn + u)/N]\ includes the photon num¬ 
ber operator n and accounts for a realistic detector with 
quantum efficiency r] and dark counts v. In this expres¬ 
sion, the exponential operator corresponds to the pro¬ 
jector of the vacuum operator, and as expectation value 
yields the probability of zero clicks [33]. 

In order to quantify the binomial character of the click¬ 
counting statistics, we employ the parameter 

Qb = h with (1) 

{ > 0 for super-binomial click statistics, 

= 0 for binomial click statistics, 

> 0 for sub-binomial click statistics, 

where (k) = J2fc=o ^Ck and ((A k) 2 ) = ( k 2 ) — ( k ) 2 rep¬ 
resent the average number of clicks and the variance 
thereof, respectively [34] • For an ideal coherent state, 
the mean and the variance of the click-counting statistics 
are analytically calculated as (k) = Np and ((A k) 2 ) = 
TVp(l— p), respectively, withp = 1 — exp[— (? 7 |<a| 2 + zA/TV]. 
When substituting these two expressions into Eq. ^ we 
find that Qb — 0 regardless of the quantum efficiency and 
the dark count rate. Accordingly, it is a sufficient crite¬ 
rion to differentiate between classical and non-classical 
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light: While super-binomial click statistics such as these 
of realistic laser light are characterized by Qb > 0, gen¬ 
uine quantum entities, such as Fock states, necessarily 
features Qb < 0. Note that in order for this parameter 
to be meaningful, the ensemble has to include N > 2 on- 
off detectors, otherwise measurements of any input state 
would produce Qb = 0 [34 , S3 • 

In a first set of experiments, we measured the absolute 
number of click coincidences M & and the relative frequen¬ 
cies = Mk/M for attenuated laser light, where k 

again represents the number of clicks within a time win¬ 
dow of 10 ns and M = M & is the total number of 

time windows. Further details on our setup are given in 
the Appendix [Aj From these measurements, we extract a 
positive value Qb = (1.712±0.026)*10 -2 , confirming that 
the click statistics is indeed super-binomial as expected 
for a classical coherent state. Beyond the Qb parame¬ 
ter, which relies solely on second-order correlations, one 
can also employ the higher-order correlations contained 
within the matrix of moments to identify non-classical 
behavior (see Appendix H- 

In our second set of experiments, we used our de¬ 
vice to characterize the fidelity of a heralded single¬ 
photon source based on spontaneous parametric down 
conversion. As in the previous case, the input 
state l^in) = af |0) is spatially distributed and thus 
transformed according to the expression |Vy>ut) = 
(— a\+\a\+a\+\a\— aj+iag— a\— idg)|0)/\/8, where aj 
represents the bosonic creation operators of the / th 
waveguide mode. The high fidelity of our device is con¬ 
firmed by the homogeneity of the single photon num¬ 
ber output distribution with an average of 12.5 =b 0.6 
(see Fig. [lj bottom). Measurement data clearly demon¬ 
strates the non-classicality of the input state, with Qb — 
-(2.002 ± 0.049) • 10“ 2 < 0. 

Finally, we analyzed the dependence of the fidelity of 
our device for different source brightness levels. To this 
end, we determined the Qb parameter for classical laser 
light with different attenuation ratios, as well as for Fock 
states at different count rates. Figure [3] illustrates that 
our scheme allows for a clear distinction between clas¬ 
sical and non-classical behavior. Note that, as the in¬ 
cident intensity decreases, the chance of multiple pho¬ 
tons entering the device in any given time slot gradu¬ 
ally converges to zero, as does the corresponding Qb pa¬ 
rameter. Whereas the sign of Qb remains well defined 
throughout this process, the uncertainty determined by 
the experimental conditions [see Appendix [A] eventually 
exceeds the absolute value. Nevertheless, standard quan¬ 
tum sources routinely feature count rates well above this 
limit, placing them firmly within the window of confi¬ 
dence of our characterization scheme. 

In conclusion, we have introduced a new paradigm for 
integrated photon-number resolved measurements based 
on distributed sensing with multiplexed arrays of con¬ 
ventional on-off detectors. In this divide-and-conquer 
approach, the limiting factor of detector dead time is 
overcome by transforming the incident fields into ex- 
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FIG. 3: Experimentally obtained Qb for two different click 
statistics versus a variation of the average click number (k). 
The inset shows a semilogarithmic section for the smallest 
achieved (k) s. An attenuated laser light (blue curve) is always 
accompanied by a positive Qb (super-binomial), whereas our 
single photon source (green curve) constantly shows a neg¬ 
ative Qb (sub-binomial). In both cases, a coherent state 
(binomial photon number distribution) is covered within the 
error-bars at very low count rates. 

tended uniform distributions. As such, coincidences in 
the same measurement channel are reliably suppressed 
even for highly multi-photon input states. Measuring 
the click statistics of the detector ensemble therefore pro¬ 
vides the means to determine the actual photon-counting 
statistics. Consequently, our technique can in principle 
be scaled to allow for arbitrarily high numbers of inci¬ 
dent photons, irrespective of the dead time of the indi¬ 
vidual detectors used. On a higher level, similar mul¬ 
tiplexing schemes could even be employed to boost the 
performance of PNR detectors relying on other mech¬ 
anisms of registration. Our experiments pave the way 
towards genuine integrated photon-number-resolving de¬ 
tectors for advanced on-chip photonic quantum networks. 

Appendix A: Methods 
1. Fabrication &; characterization 

The multiplexing waveguide network was fabricated 
by means of the direct femtosecond laser inscription in 
fused silica glass [29, 30 . It is designed to match the 
technical standards of the attached V-groove fiber ar¬ 
ray with 127 /im pitch which collects the photon outcome 
and feeds it into single-photon click-detectors. In order 
to inject weak laser light, different neutral density (ND) 
filters were placed into the beam of a laser diode emit¬ 
ting at 808 nm. The attenuated light was coupled into 
a single mode (SM) polarization maintaining (PM) fiber 
attached to the input wave guide of the multiplexer. For 
the other set of experiments, single photons of 815 nm 
were generated by spontaneous parametric down conver¬ 
sion (SPDC) in a BiBsOe crystal and coupled into two 
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SM and PM fibers. One of the fibers was directly con¬ 
nected to a single-photon detector to provide a herald 
for the other twin photon, which was delivered directly 
to the device’s injection site. By using the second pho¬ 
ton as trigger, we can suppress any noticeable influence 
of dark counts. For reasons of experimental convenience, 
the overall photon flux was again controlled by an appro¬ 
priate choice of ND filters. 


2. Measurement methods &; errors 


The photon clicks by our APDs were collected by a 
time tagging card capable to handle up to 16 detectors 
at the same time. The coincidence time window, wherein 
two or more clicks are interpreted as a joint event, was set 
to At = 10 ns for all measurements. In turn, the num¬ 
ber of measurements for attenuated laser light was deter¬ 
mined by T/At where T is the overall recording time. In 
contrast, the number of quantum measurements was de¬ 
termined by the photon flux of the trigger photons, and 
therefore remains independent of the coincidence time 
window. In both cases, the overall measurement time 
was chosen in a way to collect at least 100 million non¬ 
zero click events, ranging from 44 seconds up to more 
than 60 hours. 

The accuracy of the calculated Qb parameter is chiefly 
dependent on the number of measurements. In an ideal 
set of data, the frequencies Ck would be distributed ac¬ 
cording to the graphs shown in Fig.[2j If a measurement 
is terminated before the highest number of joint clicks 
yields a non-zero quantity, it cannot perfectly match the 
(sub- or super-) binomial distribution. Consequently, the 
extracted Qb parameter is subject to a systematic un¬ 
certainty. A second source of errors results from imper¬ 
fections in the splitting ratio of the multiplexing device, 
corresponding to small deviations from the ideal case of 
homogeneously distributed outputs (see Fig. [I]). This 
also influences the coincidence click statistics and, hence, 
the error for the Qb parameter. In turn, the resulting 
overall range of error defines a certain minimum count 
rate. 

We would like to emphasize that our method is capable 
of reliably identifying the signature of the input state, 
whether classical or quantum (see Fig. [3]), without the 
need for any corrections or post processing of raw data. 


Appendix B: High-order correlations 


In order to analyze high order correlations of the light 
fields under study, we use the approach recently intro¬ 
duced in [33] . Along these lines, we compute the matrix 
of moments MoM directly from the experimental click¬ 
counting statistics C^ xp . This is done using the expres¬ 


sion MoM = ((:d m+n :)) m?n . Here, 

(:t r m+n :) (Bl) 

= {N_Jm+n\y. ... (fc-[m+n]+l)Cf p , 

k =0 

with the indices m,n ranging from 0 to [N/2 J, where 
|_.. .J represents the floor function. It can be shown that 
the second principal minors of the MoM can be ex¬ 
pressed in terms of the Qb parameter, and as a result 
for classical states the associated MoM is expected to 
be non-negative. The positive semi-definiteness of the 
MoM can be established through the analysis of corre¬ 
sponding eigenvectors and eigenvalues. In other words, 
the MoM will be positive definite if all its eigenvalues 
are positive. Therefore, if we find at least one negative 
eigenvalue, we can unequivocally conclude that the light 
field being characterized is non-classical. This condition 
is mathematically described by 

LAT/2J 

0 > pMoMf = (:/t/ : ) = ]T 

m,m '=0 

(B2) 


where we have defined the operators / = /o + /i7r + ... + 
/ljv/ 2 J^ LJV/2J and /= (/o,..., /[iv/ 2 j) T represents the 
eigenvectors of MoM. For a classical coherent source, 
the moments can be calculated analytically using the ex¬ 
pression (:7r m :) = p m = (1—exp[—(? 7 |a| 2 + z/)/7V]) m . This 
indicates that the of an ideal coherent state exhibits one 
positive eigenvalue given by (1 — _p 2 L iV / 2 J)/(1 — p 2 ), and 
a total of [N/2\ zero eigenvalues for N > 2. 

Here, we consider a l-to-8 multiplexer device, which 
corresponds to a 5 x 5 MoM: 


MoM = 


p°:> 

W 4; > 


(:ir s :)J 


(B3) 


Using the measured C)l xp , we numerically obtain the 
eigenvalues and eigenvectors. The results of this analysis 
for experimental data obtained over 83.2 s (corresponding 
to windows) are shown in Tables [T] and [TTj for the classi¬ 
cal and quantum states, respectively, whereas Figure [3] 
shows a graphical representation of these quantities. 

The significance for these experimentally obtained 
eigenvalues can be quantified using the expression E = 

°-(0/V=>) ’ where is the standard error of the 

mean (:p f:). Note that the two negative eigenvalues, 
which unexpectedly occur in the classical case, can be 
considered zero since their absolute values fall within 
the error margin. The values of E for all measure¬ 
ments are shown in the fourth column of Table [IJ il¬ 
lustrating that all eigenvalues are non-negative within 
the range of their significance. A similar validity check 
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0.0976064 

4 

-8.59786 • 10“ 16 ±3 • 10“ 15 

8.87048-10 -15 

0.0969266 

5 

-2.47982 • 10“ 22 ±3 • 10“ 22 

6.09801 • 10“ 22 

0.406661 


TABLE I: Correlation measurements for classical coherent light. As expected from a classical source no one of the eigenvalues 
is significantly negative. 


Eigenvalues 
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841.894 
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4.47567 • 10“ 14 

0.845747 

4 

2.78037- 10 -21 ±3- 10“ 15 

3.29078 • 10“ 21 

0.844898 

5 

0±0 

0 

Indetermined 


TABLE II: Correlation measurements for the single-photon Fock states. 


was also applied for our measurements of single pho¬ 
ton Fock-states (Table 0 - As in the previous case, we 
calculate the significance of the deviation from the clas¬ 
sical bound, i.e. zero. This is done using the second 
eigenvector / = ( 0.00384503, -0.999993, -0.000251841, 
—5.430982775 • 10 -8 , 0) T which determines the operator 
test. The obtained value of E = 841.894 certifies that 
the state under consideration is indeed non-classical. 
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